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We analyze the magnitude-redshift data of type la supernovae included in the Union and Union2 
compilations in the framework of an anisotropic Bianchi type I cosmological model and in the 
presence of a dark energy fluid with anisotropic equation of state. We find that the amount of 
deviation from isotropy of the equation of state of dark energy, the skewness 5, and the present level 
of anisotropy of the large-scale geometry of the Universe, the actual shear So, are constrained in 
the ranges -0.16 <5< 0.12 and -0.012 < So < 0.012 (Ict C.L.) by Union2 data. Supernova data 
are then compatible with a standard isotropic universe (J = So =0), but a large level of anisotropy, 
both in the geometry of the Universe and in the equation of state of dark energy, is allowed. 

PACS numbers: 98.80.Es, 98.80.Jk 
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I. INTRODUCTION 



Probing the large-scale geometry of the Universe at 
cosmological scales is, undoubtedly, one of the most 
outstanding issues in modern cosmology. The stan- 
dard assumptions of homogeneity and isotropy (namely, 
the Cosmological Principle can now be tested via 
new and very accurate data coming from the study of 
cosmic microwave background (CMB) radiation, espe- 
cially from the Wilkinson Microwave Anisotropy Probe 
(WMAP) [2], and data from type la supernovae, such as 
those collected in the so-called Union [3] and Union2 [4] 
compilations. Indeed, concerning the tests of isotropy, 
an anisotropic model of Universe, known as "ellipsoidal 
universe" |5l4ll|. may be even favored by an observed 
anomalous feature of the CMB power spectrum -the lack 
of power on large angular scales- while being consistent 
with other cosmological data. 

It is not excluded that such an anisotropic model of 
the Universe could even account for three other large- 
scale "anomalies" of the isotropic standard cosmological 
model (for a brief but pointed discussion see Ref. [12|): 
the detection of large-scale velocity flows significati vely 
larger than those predicted in standard cosmology |13| . 
a statistically significant alignment and planarity of the 
CMB quadrupole and octupole modes [lj|, and the ob- 
servation of large-scale alignment in quasar polarization 
15|. It should be stressed, however, that the 



vectors 

above large-angle anomalies are still subject to an in- 
tense debate, since they could be indeed related to some 
common systematic. 

In this paper, we present an analysis of the large-scale 
isotropy assumption by means of magnitude-redshift data 
of type la supernovae (SNe). In particular, we use data 
from both Union and Union2 compilations, consisting of 



307 and 557 type la SNe respectively, to set constraints 
on the parameters of an anisotropic model of the Universe 
(for earlier work on the possibility to test the Cosmolog- 
ical Principle with SN data, see Ref. [iMi, while for 
recent works, see Ref. 2^ 21 1 and references therein). 



We assume an anisotropic Bianchi type I cosmologi- 
cal model [i^], characterized by a cosmic shear S in the 
presence of a dark energy fluid with anisotropic equa- 
tion of state, characterized by a skewness 5. This fluid, 
firstly studied by Barrow in Ref. 23], could be produced 
by the dynamics of a cosmic vector field, as shown by 
Koivisto and Mota in Ref. [l^. Other effects could give 
rise to an ellipsoidal universe, such as a large-scale cosmic 
magnetic field , or a dark energy fluid having a non- 
vanishing velocity with respect to the CMB frame p^ . 
For an incomplete list of such mechanisms of universe 
anisotropization see, e.g., Ref. 1^ and references therein. 

Testing the dark energy anisotropic model with type 
la SNe implies that we can only constrain the anisotropy 
parameters (S and 5) at relatively recent times (i.e., at 
redshift z < 1.6), their earlier evolution being largely 
unconstrained. 

We find no evidence in favor of anisotropics of either 
geometric origin (E ^ 0) or dark-energy origin {S ^ 0). 
However, we can put significant upper and lower bounds 
on the deviations of E and S from zero. 

The paper is organized as follows. In the next Sec- 
tion we set up the formalism of a cosmological model 
with anisotropic fluid while, in Section III, we derive the 
magnitude-redshift relation for such a universe. In Sec- 
tion IV, we use magnitude-redshift data of SNe from the 
Union and Union2 compilations to constrain all the free 
parameters of the model, including E and S, so as to test 
the isotropy of the observable universe at z < 1.6. In 
Section V, we draw our conclusions. 
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II. ANISOTROPIC COSMOLOGICAL 
MODEL: ELLIPSOIDAL UNIVERSE 

In order to test possible anisotropies of the Universe 
we need to make assumptions beyond the standard cos- 
mological model (the Lambda cold dark matter (ACDM) 
concordance model which is isotropic. 

For the sake of simplicity, in this paper we consider 
an anisotropic Bianchi I type cosmological model with 
the highest (planar) symmetry in the spatial sections of 
the spacetime geometry. Named as "ellipsoidal universe" 
in (Hl-lllt. it has the attractive feature of accounting for 
the observed lack of power of the cosmic microwave back- 
ground anisotropy at large scales. 

In this model, the most general plane-symmetric line 
element is 1251: 



'-{t){dx'^ +dy'^)-b^{t)dz^ , (1) 



where a and b are the scale factors (or expansion param- 
eters), which we normalize as a(to) — 6(io) = 1 at the 
present time Iq. ^ The mean Hubble parameter H is 
defined as 



H 



A 

A 



(2) 



A = {a^b)^^^ being the mean expansion parameter, while 
the Hubble parameter in the symmetry plane is 



a 



(3) 



Here and in the following a dot denotes differentiation 
with respect to the cosmic time. The cosmic shear S is 
defined as: 

H ^ ' 

In an ellipsoidal universe, the most general energy- 
momentum tensor compatible with the metric in Eq. ([T]) 
is of the form 



diag(p, -P||,-P||,-Pi) , 



(5) 



where p\\ and pj_ are "longitudinal" and "normal" pres- 
sures. In Fricdmann universes, p\\ — p±, reflecting the 
isotropy of the metric. Conversely, anisotropic universes 
can support fluids with p|| ^ which we will refer to 
as anisotropic fluids. 

Given the above energy-momentum tensor, Einstein's 
equations read 

3(1 - S2)ij2 = SttGp , (6) 
3(1 - S + E2)i/2 + £[{2- S)i7] = -8^Gp|| , (7) 

3(l + I])2i/2 + 2^[(l + I])i7] = -SttGp^ , (8) 



where G is the Newton constant. 

By linearly combining Einstein's equations we get 

-(H^) + 3H^j:^—ip\\-p^), (9) 

which shows that asymmetric pressures can act as a 
source of the shear, and 

P + 3h(^P+ ^^11 ^ ) + 2g(p|| - p^)S] = , (10) 

which represents the time component of the energy- 
momentum conservation law, T^^.^ — 0. 

We assume that the anisotropic fluid defined in Eq. ([S]) 
is indeed made up by an isotropic, pressureless dark mat- 
ter (DM) component and an anisotropic, dark energy 
(DE) component, with equations of state 

p\l = W|| pde , P± = wj_ pde , (11) 

where pde is the dark energy density. For the sake of 
simplicity, we assume constant and w_l coefficients. 
In terms of and w±, we define a mean coefficient w 

as 

2w\\ + w± 



and a skewness 6 as 



S = w\\ — w± 



(12) 



(13) 



We also assume noninteracting DM and DE fluids, so that 
the isotropic and anisotropic components are separately 
conserved. Writing 

P = Pm + PUE , (14) 
where pm is the dark matter density, Eq. ((TU)) gives 

p„ = p^^jA-' , (15) 
PDE + [3(1 + w) + 2ST,] HpuE - , (16) 

where the superscript "0" denotes quantities evaluated 
at the present time to- 

We finally recast Einstein's equations in terms of di- 
mensionless quantities as 



HA 
A 



d{H^) 



dA 



3i/^E = ^r^DEPDE , 



(17) 



with 



dA 



H = 



[3(1 + w) + 2<5E] PDE = , (18) 



ilmA 3 + I^DEPDE 



1-S2 



(19) 



In Eqs. P?)) -(|19 p . each barred quantity is normalized to 
its actual at to, and the usual energy density parameters 
are introduced, 



^ In the standard cosmological model, a = fe at all times. 



(0) (0) 

-(Oj ' ^^DE = — , (20) 
Per Per 



3 



the critical density being at to 



which is connected to the shear via Eq. ^ as: 



P. 



(0) ^ 3gg 
~ SttG 



(21) 



It is worth noting that, given fi™ and Sq, f^DE is derived 
from Eq. as 

f^DE = 1 — i^m — Sq . (22) 

For given parameters {V,„i,w,S,'Eo}, Eqs. PT)) "([T ^ are 
solved numerically. The resulting solution, once the value 
of Ho is fixed, describes the evolution of the ellipsoidal 
universe. 



III. DISTANCE MODULUS IN ELLIPSOIDAL 
UNIVERSE 

Within the anisotropic metric ([T]), the luminosity dis- 
tance of a source at redshift z, seen along the direction 
p i\p\ = 1), is given by 



dL{z,p) = (1 + z) 



dt 



where ai = a2 = a and 03 = 5 and [3, [2^, H^l 



(23) 



(24) 



The above result applies to the case where the axis 
of symmetry is identified with the z-axis. To account 
for a symmetry axis directed along a generic direction 
(6, 1) = {hA, I a) in galactic coordinates [b and I being the 
galactic latitude and longitude, respectively) , we perform 
rotations of the coordinate system along z and x, 



n = n47r/2 - bA) 7^^(7r/2 + Ia) 



(25) 



where the arguments are the rotation angles. In the 
galactic coordinate system, the direction cosines of the 
symmetry axis are 

fiA — (cos6yi cos Z^, cosfoyi sinZ^, sin6yi) , (26) 

while the components of the generic direction p are de- 
fined by direction cosines h — TZ^^p, namely, 

n = (cosfocosZ, cosfosinZ, sinfo) . (27) 

The angle 9 between h and fiA is defined by 

cos6' = n ■ fiA ■ (28) 

It is useful to introduce the "eccentricity" e as 

&2 



6^ = 1 



(29) 



1 — exp 



6 ^ I](x) 

J A X 



(30) 



In galactic coordinates system, the redshift and dis- 
tance modulus read then 



1 1 - sin^ 



A 



(l-e2) 



2\l/3 



(31) 



and 



dL{z,e) 



\^z 



dA 



Ja(z) ^^i^ (l-e2cos2 



.1/2 



(32) 



where H and e are taken as function of A, and A{z) is 
the solution of Eq. (|5T|) . 

Finally, we introduce the usual distance modulus n, as 



^ = 5 log 



i°(lMpc) 



25 . 



(33) 



It depends on two different sets of parameters, fj, — 
where i = {z,b,l} characterizes the position of 
the supernova, while j — {Hq, ftm, Sq, w, S, bA, Ia} char- 
acterizes the cosmological model. 

In the next section, we use the magnitude- redshift SN 
data to constrain the free parameters of the the ellipsoidal 
universe model with anisotropic dark energy. 



IV. TYPE IA SUPERNOVAE: 
CONSTRAINING THE LEVEL OF COSMIC 
ANISOTROPY 



In a given cosmological model, the analysis of the SN 
lightcurves provides, for each supernova, the coordinates 
(z,/i) which form the so-called Hubble diagram. 

In the following, we discuss the Hubble diagram for the 
Union compilation Q , which consists of 307 supernovae, 
and the recently published Union2 compilation con- 
sisting of 557 supernovae. 

For our anisotropic cosmological model, the distance 
modulus depends (besides the cosmological parameters) 
on the redshift and angular position of the source. Hence, 
in order to fit experimental data, we need for each super- 
nova its angular position. Figure 1 shows the angular 
position (6, 1) of SNe in the Union and Union2 compila- 
tions, b and I being the galactic latitude and longitude, 
respectivel y. Data on angular positions are taken from 
Ref. [1, [isl - liol . The non- uniform distribution of the data 
(e.g. the absence of SNe around the galactic plane at 
b = 0) reflects observational selection cuts which, how- 
ever, do not affect our subsequent fit procedure. 
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FIG. 1: Angular position of type la supernovae in the Union (left panel) and Union2 (right panel) compilations, h and I are 
the galactic latitude and galactic longitude of supernovae, respectively (in degrees). 



IV. a Union Compilation 

From the observed SN lightcurve, one can infer the 
distance modulus /Lt^ as the difference between the peak 
bolometric apparent magnitude m^^^ and the absolute 
magnitude M of supernova. The derived distance mod- 
ulus, as well known, is affected by the so-called stretch 
and color corrections, s and c respectively. The amount 
of such corrections is adjusted by means of free parame- 
ters a and /3 as [sj 



(34) 



with M, a and (3 to be determined from fits to data. The 
experimental values rw^^^, s, and c for each supernova 
are taken from Ref. [3]. 

In order to constrain the free parameters of our 
anisotropic cosmological model, we compare the exper- 
imental distance modulus IjLb{M, a, jS) with the theoreti- 
cal expectation ii{Hq, Q.^, Eqj Wi 5, bA, I a) by means of a 
least-square fit 



X 



E 

SNo 



'AM, 



(35) 



which we now discuss in detail. 

Theoretical systematics. The uncertainty cr^ of the 
theoretical distance modulus /i contains two independent 
contributions, 



lightcurve 



(36) 



The first one, uiightcurvc, is associated to the fitting 
procedure of SN lightcurves and can be expressed as 



Si 7 C-iCjCij 



where Ci = {l,a, — /3} and 



lightcurve 

Cij is the covariance matrix, whose diagonal parameters 
are the uncertainties on peak bolometric apparent mag- 
nitude, CTto^, stretch, ct^, and color ct^. Since the co- 
variance matrix for each supernova, obtained from the 
lightcurve fitting procedure in Ref. 9], is not publicly 
available, we assume uncorrelated errors on m^, s, and 
c, but a nonzero stretch-color correlation from Ref. j4l| , 
Psc — —0.2. Accordingly, we get 



lightcurve 



a'ai 4- 13' ai + 2p,, a ^ a^ac ■ (37) 



The uncertainty (Tastro is made up of three different un- 
correlated contributions of astrophysical origin Q : 



o-ie 



(38) 



The first (av-pec) is associated to the peculiar velocities 
Wpec of the host galaxies, and is estimated by generalizing 
the analysis of Ref. 42[ (valid in a standard isotropic 
universe) to the case of an ellipsoidal universe. We find: 



where 



with 



V — pec f {.^1 ^) ^pCC 5 



f{z,e) = i 



HdL 



(39) 



(40) 



(1 - e^)V^ (1 - sin"g)i/^ (1 - cos^g)'!/^ 
(1 - e2 sm'^9){l -t- S) - 3E cos^9 ' 

(41) 

Equation (15^]) reduces to the result of Ref. [42| in the 
limit of vanishing eccentricity e. As a typical peculiar 
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FIG. 2; \J Ax^ for 9 of the 11 parameters of the anisotropic cosmological model using the Union set. The two probabihty 
distributions for I a and hA are flat (not shown). 



velocity, we take the value fpoc — 300km/s. 

The second term (aicns), is associated to gravitational 

lensing and can be parameterized as j3j] 



0.093 z 



(42) 



Finally, the third term (ccxt) is associated to Galactic 
extinction corrections, and can be estimated as: ^ 



fext — 



0.013 , z < 0.2 
0, z>0.2 



(43) 



In Ref . 01 , the extinction law of Cardelli et al. [43j , together with 
the E{B — V) values derived from the sky map of Schlegel et 



al. [44j . were used for parameterizing Galactic extinction, 
we just use that law with the average Galactic E{B — V). 



Here, 



Experimental systematics. In order to take into account 
systematic errors on fis, we follow the analysis of Ref. [3] . 
Firstly, we introduce a sample-dependent uncertainty 
CTsys related to an unknown intrinsic dispersion of the 
supernova magnitudes, in each of the 13 data subsets in- 
cluded in the Union compilation. For completeness, these 
uncertainties are reported in Table I. 

Secondly, we add two nuisance parameters, AMi and 
AM2, shifting fiB as 



AMi 
AMi 



AM, 



z < 0.2 , 
z > 0.2 . 



(44) 



The nuisance parameter AM2 is z-dependent, since it 
accounts for a possible evolution of supernovae with red- 
shift. The value z = 0.2 discriminates, conventionally, 
between "low" and "high" redshifts in this context [3]. 
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FIG. 3: la, 2a, and 3a confidence level contours in the planes (Eq, S), {w, 5), {flm, S), (w, Eo), (f2m, Eo), and (fim, w). Results 
obtained by using the Union set are shown as gray filled contours, while those coming from Union2 are represented as empty 
contours. 



The uncertainties on the nuisance parameters, crzv\/i and 
f AM2 J are found by adding in quadrature various uncor- 
related uncertainties on the a and /3 parameters, on con- 
tamination, on hghtcurve model, on zero point, on the 
Malmquist bias, and on galactic extinction. Taking into 
account the results of Ref. 3| we use the values: 



CTAA/i = 0.040 , (JAM, = 0.034 



(45) 



The previously introduced function is then completely 
defined. 

Notice that the two parameters Hq and M enter in the 
only through a specific combination Hm, 



M 



km/s/Mpc 



-^Q-(Af+25)/5 



(46) 



Therefore, in the analysis we cannot constrain H and 
M separately, but only Hm- We expect Hm — 5.1 for 
small deviations from standard cosmology, where Hq ~ 
71km/s/Mpc and M ~ -19.3. 

Since the depends on as many as eleven parameters, 
a brute- force minimization search on a 11-dimensional 



grid is not feasible. We employ then a Markov Chain 
Monte Carlo approach, which has become standard in 
CMB analyses [3] • We use a modified version of the Cos- 
moMC (Cosmological Monte Carlo) code {isj to produce 
and analyze the likelihood chains. 



Results. Figure 2 shows the \/Ax^ dis- 

tributions for nine parameters of the model: 
w, r2„i, i?M, So, (5, a, /3, AMi, AM2. The two remaining 
parameters I a and Ba, which define the symmetry axis, 
are basically unconstrained. This means that there is 
no evidence for a specific "anisotropy axis" from the 
data (see also discussion in Section IV. c). However, for 
any chosen direction, the data can constrain the level of 
anisotropy, characterized by the parameters Eq and S. 

In Fig. 3 (see gray filled contours), we present the la, 
2(7, and 3a contours in the planes (Eq, S), {w, 5), (fim, 5), 
{w,Tio), (OmjEo), and (flm , w) . Here again we follow 
the convention of Ref. [49|, where the la, 2a, and 3a 
joint regions are defined so that the projections onto each 
parameter give, respectively, the la, 2a, and 3a intervals 
for that particular parameter. The at the best-fit 
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TABLE I: Intrinsic dispersion of the supernova magnitudes, 
CTsys, for each subset of the Union compilation. The second 
column indicates the number of supernova in each subset. 
Data from Ref. H. 



Subset of Union 


SNe 


*^sys 


Ref. 


Hamuy et al. 


17 


0.14 


[28] 


Krisciunas et al. 


6 


0.05 


[39] 


Riess et al. 


11 


0.16 


[45] 


Jha et al. 


15 


0.26 


[29] 


Kowalski et al. 


8 


0.00 


[3] 


Riess et al. + HST 


12 


0.28 


[45] 


Perlmutter et al. 


29 


0.33 


[46] 


Tonry et al. 


6 


0.06 




B arris et al. 


21 


0.23 


[32] 


Knop et al. 


11 


0.10 


[47] 


Riess et al. 


27 


0.16 


m 


Astier et al. 


71 


0.12 


[37] 


Miknaitis et al. 


73 


0.18 


[38] 



point is 255.2. 

An isotropic universe ((S = Eq = 0) is consistent with 
data and present uncertainties and statistics only allovif 
us to put upper limits on the level of cosmic anisotropy. 

The most correlated variables are w and flm (bottom- 
right panel), and S and Eq (top- left panel). 

The anti-correlation between w and Qm is also present 
in the standard analysis of SN data [3] and it comes from 
the peculiar dependence of luminosity distance on w and 
rim. If one combined SN data with other cosmological in- 
formation, such as CMB or Large Scale Structures anal- 
yses 1], that degeneracy could be broken so to better 
constrain S and Eg. However, as it is clear from Fig. 3 
[see in particular the {w,S) and (r2„j,(5), and the (w,Eo) 
and (r2,„,Eo) panels], the bounds on S and Eq depend 
only weakly on w and Q^n- 

That S and Eg should be correlated was to be expected 
since the deviation from isotropy of the equation of state 
of dark energy is the source of the anisotropization of 
the large-scale geometry of the Universe. Indeed, in the 
limit of small anisotropics, i.e. \S\ <^ 1 and |E| ^ 1, the 
correlation between 6 and Eg can be understood as fol- 
lows. To the first order in the anisotropics, the solutions 
of Eqs. (HZl) and dH]) are 



E(A) = 



Eq + (^ - ^o) S 



and 



-3(l+-!«) 



Pde = A 

respectively, with the Hubble parameter being 



(47) 



(48) 




0.4 0.6 0.8 1.0 1.2 1.4 1.6 



FIG. 4: Upper panel. Hubble diagram for the 307 supernovae 
in the Union compilation, for different cosmological models: 
best-fit cosmology (Eo, S) ~ (—0.012, —0.013) (red continuous 
line), (So, (5) — (0.2, 1) (magenta dash-dotted line), (Eo,(5) = 
(0.2,0) (green dashed line), (Eo,5) = (0.2,-1) (blue dotted 
line). For graphical clarity, we take the SN angular positions 
and the direction of the axis of symmetry fixed to the value 
{b,l) — (&a,'a) ~ (0,0). The remaining parameters are fixed 
to their best-fit values. Lower panel. Residuals (distance 
modulus minus distance modulus for the best-fit cosmology) 
for the same models in the upper panel. 



In Eq. ([47]) we have introduced the function 

/■^ dx 
Jo 

and the quantity 



H{x) 



(50) 



Eo 




where B{x; a, b) is the Incomplete Euler Beta Function 
and the last equality is valid for w < 0. 

At early times {A ^ 1), the anisotropic source as- 
sociated to dark energy is inefficient to generate an 
anisotropization of the universe since the dark energy 
density is completely negligible with respect to that of 
the isotropic dark matter component. This condition of 
isotropization at early times, namely 



(49) 



lim E(A) 







(52) 





-0.10 -0.05 0.00 0.05 0.10 -1.0 -0.5 0.0 0.5 1.0 -3 




FIG. 5: \J Ax^ for 4 of the 6 parameters of the anisotropic cosmological model using the Union2 set. The probabihty distribu- 
tions for I A and 6a are flat (not shown). 



when applied to Eq. (|T7| gives 



IV. b Union2 Compilation 



So 



(53) 



Observing that is a positive-defined quantity which 
depends only on w and Sim, we find that & and Eq are, 
indeed, positively correlated. Eq is of order of 10~^ for 
a fiducial cosmology (Sl„i,w) — (0.3,-1). So, we expect 
(5 ^ lOSoi and this is confirmed by our data analysis (see 
the top- left panel in Fig. 3 and, below. Table II). 

Finally, we observe that the a and /3 best-fit values 
are in good agreement with those found in Ref. Q for an 
isotropic universe, and that the probability distribution 
for the nuisance parameters AMi and AM2 are centered 
around zero, indicating that no systematic shift of the 
distance modulus is needed to fit data. Also, as expected, 
we find Hm — 5.1. 

For completeness, we present in Fig. 4 (upper panel) 
the Hubble diagram for different realizations of our cos- 
mological model. For graphical clarity, we take the SN 
angular positions and the direction of the axis of sym- 
metry fixed to the value (6, V) = (6^, ^a) = (0, 0). In the 
lower panel we show the corresponding distance modulus 
fit residual (distance modulus minus the best-fit distance 
modulus) as a function of the redshift. The red contin- 
uous line corresponds to our best-fit cosmology, practi- 
cally indistinguishable from the one corresponding to the 
standard isotropic universe (not shown) . For a relatively 
large value Eo = 0.2 (that we choose for illustrative pur- 
poses) and for (5 = 0, the distance modulus curve moves 
downwards (green dashed line). By taking (5 = 1, the 
distance modulus curve moves upwards (magenta dash- 
dotted line), while the opposite happens taking 6 — —\ 
(blue dotted line). A similar behavior is found for nega- 
tive value of the shear, according to the positive correla- 
tion found between Sq and b. 



The analysis of SN lightcurves in the Union2 compila- 
tion gives the distance modulus [Ib as 3 

^is = mg^'^ - M olxx - /3'c , (54) 



where rn^^'^ 



and c are, as before, the peak bolometric 
apparent magnitude and the color correction, while xi 
derives from the fit to the SN lightcurves with the so- 
called SALT2 (Spectral Adaptive Light curve Template) 
fitter a' and /3' are, instead, nuisance parameters 

to be determined, simultaneously with the cosmological 
parameters and the absolute magnitude M, from fits to 
data. 

statistics. As in the case of Union compilation, the 
covariance matrices resulting from the lightcurve fit are 
not publicly available, so we do not have any information 
on the correlation between the errors on m^^^, xj, and 
c. Consequently, we follow the analysis of Ref. [48| and 
introduce the statistics as 



X = mm 

7 



(55) 

where the double sum runs over the 557 SNe, /i^'^^ is the 
experimental value of the distance modulus of the i^^ su- 
pernova, Hi the corresponding theoretical value given by 
Eq. ([33]) . afj is the covariance matrix (containing both 
statistical and systematic errors), and 7 an unknown nor- 
malization parameter 5l| . The values of hI^^ and afj in 
the our analysis are taken from Ref. 

The minimum in Eq. (|55p can be evaluated analytically 
and gives 



X 



where we have introduced the matrix 

„-2 -2 

l^kl "ik "ij 



M,,. 



(56) 



(57) 



kl 
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FIG. 6; 1(7, 2(7, and 3(7 confidence level contours in the planes (Eo,5), {w,S), {Qrn,5), (w, Eo), (r2m,Eo), and {Qrn,w). Gray 
filled contours come from the Union data analysis detailed in Section IV. a, while empty contours refer to the Union2 data 
analysis of just the 293 SNe contained in both compilations. 



and the "reduced" distance modulus 

il = 5 logio dL + 25 , (58) 

with the "reduced" luminosity distance 

dL = HodL . (59) 

It is worth noticing that (Il is independent on the 
Hubble parameter Hq, so both p, and the i^i 
Eq. (|56l) depend only on the cosmological parameters 
(fim, So, w, (5, 6a, Z^). Since the depends on as many 
as six parameters, we employ a Markov Chain Monte 
Carlo approach. 

Results. Figure 3 (empty contours) shows the la, 2a, 
and 3(7 contours in the planes (I]o,<5), iw,5), (fim,5), 
(uIjSq), (51m, So), and {flm,,w)- Figure 5, instead, shows 
the Ax^ distributions for So, and flm, the two 

distributions for I a and bA being flat (not shown). Fi- 
nally, Table II shows the best-fit values and the 1 and 2a 
confidence level intervals for the same parameters. The 
X^ at the best-fit point is 528.9. 



TABLE II: Best-fit (BF) values, and the la and 2a confidence 
level intervals derived from the Union2 data analysis, for the 
cosmologically relevant parameters Eo, 5, w, and flm- 



Eo 5 W Qrn 

BF -0.004 -0.05 -1.32 0.37 

1(7 [-0.012,0.012] [-0.16,0.12] [-1.66,-1.04] [0.29,0.41] 
2(7 [-0.024,0.046] [-0.29,0.246] [-1.96,-0.79] [0.22,0.46] 



The Union2 data analysis gives more stringent bounds 
on the anisotropy parameters with respect to the Union 
analysis. In fact, from the top-left panel of Fig 3, it can 
be seen that the uncertainties on Eo and 6 are reduced by 
about one half. Moreover, going from Union to Union2 
analysis, the uncertainty on w is practically unchanged, 
while slightly higher values of flm are preferred. 

However, both the anti-correlation between w and flm 
and the correlation between S and Eo still persist, as well 
as the weak dependence of S and Eo on w and flm- 
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0.50 
0.76 
1.00 




to 

T3 
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0.00 
0.26 
0.50 
0.76 
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0.8 
Z 



FIG. 7: Left panels. Upper panel: Distance modulus as a function of the redshift for different values of cos^S, with 6 being 
the angle between the anisotropy axis and the generic direction of a SN, for Eo = 0.01. The remaining parameters are fixed 
to {Ho, 5,Qrn,w) — (71 km/s/Mpc, 10Eo,0.3, —1). Lower panel: Residuals (distance modulus minus distance modulus for the 
case cos^9 = 0) for the same models in the upper panel. Right panels. The same as the left panels but with Eo — 0.1. 



In order to check if the reduction of Sq and S uncer- 
tainties is due to the increase of statistics in the Union2 
compilation with respect to Union (557 compared to 307 
SNe) or to the the different method of analysis of SN 
data, we repeated the Union2(-like) analysis keeping only 
those 293 SNe which appear in both compilations. The 
results are shown in Fig. 6, superimposed to the previous 
Union fit. The the best-fit point is 263.2. 

A part from slightly higher values of flm, and S pre- 
ferred by the 293 SN Union2 fit, the two different anal- 
yses are in very good agreement. This fact highlights 
the robustness of our bounds on the cosmic anisotropy 
obtained with type la supernovae. 

As a final check of correctness of our analysis in de- 
termining the amount of anisotropics, we verified that 
our results do not change if we exclude from the anal- 
ysis supernovae at low redshifts, which could eventually 
give spurious values of anisotropics on far larger scales. 
In particular, we cut first SNe with z < 0.02 and then 
SNe with z < 0.03, finding no appreciable changes in 
the confidence level intervals for Eq: w, and rim (not 
shown) . 



IV. c Anisotropy Axis 

We have found, in the previous Sections, no evidence 
for a preferred axis of symmetry, since the probability 
distributions for I a and Ba are flat. 



This is due to the fact that the allowed values of the 
cosmic shear Eg are relatively small, jEol S O-Ol at la 
C.L. (see Table II). 

Indeed, for small values of the shear, the distance mod- 
ulus Eq. (155)) weakly depends on both the angular posi- 
tion of SNe, h, and on the direction of the symmetry axis, 

flA- 

More precisely, these two angular positions enter in the 
expression of /i just through the particular combination 
cos^ 9 = {h ■ uaY [see Eqs. ([50]) - ([55)1 ]. and the distance 
modulus is not much sensible to the variations of cos^ 9 
when |Eo| is small. 

This can be appreciated in Fig. 7, where we show /i 
as a function of the redshift for different values of cos^fl 
in the two cases Eq = 0.01 (left panel) and Eq = 0.1 
(right panel), with 5 = lOEo and assuming the fiducial 
cosmology {Ho,ri„i,w) = (71 km/s/Mpc, 0.3, — 1). 

The variation of /i for Eq as small as 0.01 is well within 
the experimental uncertainty on the distance modulus (of 
about 0.3 magnitudes), so that a determination of the 
symmetry axis is beyond the actual experimental possi- 
bility. On the other hand, for values of the shear as large 
as 0.1 a determination of fiA is in principle possible -due 
to the significant dependence of ^ on cos^6'- but such 
large values of Eg are excluded by data. 

Therefore, we conclude that SN data are not able to 
constrain, at the same time, the anisotropy parameters 
(Eq and 5) and the preferred direction defined by the 
anisotropy itself. 
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However, it is worth noticing that the approach de- 
tailed in this paper (that provides bounds on Eq and 
5 but not on n^) can be confronted with the so-cahed 
"hemisphere comparison method" [s^] (which searches 
for a preferred axis in magnitude-redshift data without 
constraining the level of anisotropy) to give a better a 
complete understanding of the geometry of the Universe 
on large scales. 

Indeed, the hemisphere comparison method applied 
to Union2 data favors an anisotropic accelerated expan- 
sion in the direction [13] (6a, U) = (l8° tio°„ , 309° tgf ) , 

which is consistent, both qualitatively and quantitatively, 
with the evidence of a preferred direction in the sky com- 
ing from other cosmological observations (namely, large- 
scale velocity flow axis, alignment of CMB quadrupole 
and octupole axes, and quasar optical polarization align- 
ment axis, all discussed in Ref. jl2|). 



V. CONCLUSIONS 

The analysis of cosmic microwave background radia- 
tion reveals that the level of anisotropy in the large-scale 
geometry of the Universe must be below 10~^ at red- 
shifts z ^ 10'^. This constraint do not exclude that a 
much larger level of anisotropy is allowed at recent times, 
z - 0(1), as first noted in Ref. pj] and Ref. |ia]. 

In particular the authors of Ref. [l^ analyzed the 
magnitude-redshift data on type la supernovae in the 
"GOLD" data set of Riess et al. [s^] (consisting of 182 
supernovae), in a nonstandard cosmological model with 
anisotropic dark energy. They found that deviations from 
isotropy in the equation of state of dark energy were con- 
strained at the level of \S\ < few x 10^^. 

In this paper, we have extended the work [l9| by con- 
straining also the present level of cosmic anisotropy Sq, 
namely, the anisotropy in the large-scale geometry of the 
Universe, and by analyzing the magnitude-redshift data 
on type la supernovae in the Union and Union2 data sets 
of Kowalski et al. |3[ (consisting of 307 supernovae) and 
Amanullah et al. [J] (consisting of 557 supernovae), re- 
spectively. In particular, by using Union2 data, we have 
confirmed the results of [19[ about the skewness, finding 



0.16 < (5 < 0.12 (IcrC.L.) 



(60) 



and we have put first limits on the present cosmic shear 
parameter Eg, 



0.012 < So < 0.012 (iCT C.L.) 



(61) 



We conclude that a standard isotropic universe is con- 
sistent with SN data, any deviation at redshifts z < 1.6 
being constrained by the above results. 

We would like to thank E. Lisi for useful discussions 
and valuable comments. 
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